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O ■ ABSTRACT 



We find convincing observational evidence to confirm tlie optical identification of the X-ray 
burster X1746-370 located in the globular cluster NGC 6441. Chandra/ HKC-l imaging yields 
O ' a much improved X-ray position for the source, which we show to be fully consistent with our 

^ . rederived position of a UV-excess star, Ul, in the same astrometric reference frame. In addition, 

^ I the smaller Chandra X-ray error circle excludes the only other blue stars previously identified in 

the old Einstein circle. We have also obtained Hubble Space Teles cope /STIS time-resolved optical 
spectra of star Ul. Although there are no strong line features, the flux distribution demonstrates 
. Ul to be unusually bright in the blue and faint in the red, consistent with earlier WFPC2 

5-H ' photometry. More notably, the flux level of the continuum is seen to vary significantly compared 

to stars of similar brightness. Indeed, the lightcurve can plausibly be fit by a 5.73 hr period 
sinusoid, which is the period of the recurring X-ray dips seen in this source. The presence of 
modulations in both wavelengths strengthens the case for an orbital origin, and therefore deepens 
the puzzle of the unusual energy independent X-ray dips. Lastly, we note that X1746-370 remains 
the longest period confirmed X-ray burster in a globular cluster, and the only one with a period 
typical of the galactic population as a whole. 

Subject headings: globular clusters: individual (NGC 6441) — stars: neutron — X-rays: bursts — X-rays: 
stars 

1. INTRODUCTION 



1 Based on observations with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Globular clusters are expected to provide ideal 

Institute, which is operated by the Association of Universi- environments for the formation of closc binaries, 

ties for Research in Astronomy, Inc., under NASA contract ..i .1 • 1 • 1 ^ n 1 1 1 

^.T^^,^^„^^^ With tncir nigh stellar densities and much en- 

1NAS5-2d555. . . „i • • 

hanced rates 01 star interaction. This is certainly 
the case for the X-ray bright interacting systems. 
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Almost from the advent of X-ray astronomy it 
has been known that the cluster population of lu- 
minous (> lO^^ergs"^) low-mass X-ray binaries 
(LMXBs) is >100 times enhanced (per unit stel- 
lar mass) relative to the galaxy as a whole (Clark 
1975; Katz 1975). Another curious fact is that 
11 of these 13 LMXBs must contain neutron stars 
rather than black holes, as they exhibit type-I X- 
ray bursts (understood as thermonuclear runaway 
burning on the compact object's surface). More- 
over, the sensitive X-ray imaging of the Chandra 
X-ray Observatory has now revealed an equal num- 
ber of probable quiescent systems in the few clus- 
ters examined to date (Grindlay et al. 2001a,b; 
Heinke et al. 2001; Homer et al. 2001b; Pooley 
et al. 2002; Rutlcdge et al. 2002), as well as the 
existence of two persistent LMXBs in M15 (White 
& Angelini 2001). 

The study of globular cluster LMXBs benefits 
greatly from a multi-wavelength approach. X-ray 
data, in general, only probe the vicinity of the 
central source, apart from the highest inclination 
systems where material farther out can cross the 
line of sight. The periods of four cluster LMXBs 
have been determined from eclipses by the donor 
star and/or from the periodic dips in their X-ray 
flux, understood as due to obscuration by verti- 
cally extended material near the edge of the ac- 
cretion disc. However, with the identification of 
optical/UV counterparts (in all but one case re- 
quiring the resolution of HST), we can immedi- 
ately begin to estimate the linear scale of a system 
from the Lx/Lopt ratio, which has been shown to 
scale with disc area (van Paradijs & McClintock 
1994). Photometric monitoring has also proven 
effective in revealing variability on the binary pe- 
riod, whilst the optical spectra can in principle 
provide definitive corroboration of a counterpart 
and further useful diagnostics. 

The X-ray burster X1746-370, located in 
NGC6441, is one of the X-ray "dippers". From 
a continuous EXOSAT ohseivation, Parmar et al. 
(1989) first observed dips and inferred an orbital 
periodicity, which was refined to 5.73±0.15hr by 
Sansom et al. (1993) using a more extensive Ginga 
dataset. Apart from a single deep (90% fiux 
decrease) dip, which showed spectral hardening 
(Jonker et al. 2000), all the observed dips have 
been shallow (~15%) and have shown no clear 
energy dependence. The apparently energy in- 



dependent dipping is puzzling, implying that the 
obscuring material responsible for the electron 
scattering has metal abundance < 0.01 times so- 
lar, but this seems unlikely given the close to 
solar metallicity of the cluster as a whole (Djor- 
govski 1993). The alternative explanations are: 
photoionization of the material, a number of vary- 
ing spectral components conspiring together, or 
an extended X-ray source (i.e. an accretion disc 
corona). However, the most recent broad-band 
spectroscopy of Parmar et al. (1999) using Bep- 
poSAX argues against any of these possibilities. 
One might contend that the standard dipping in- 
terpretation itself could be erroneous- and it is 
certainly true that none of the X-ray period deter- 
minations has been sufficiently precise to confirm 
the recurrent period as orbital in origin based on 
its stability. Additional progress has been made 
in the optical. Using F5T/WFPC2 imaging data, 
Deutsch et al. (1998) identified a variable, UV- 
excess star (designated Ul) in the Einstein X-ray 
error circle. However, given the surprisingly large 
number of similar UV-bright stars in the cluster, 
there remained a possibility that Ul might be a 
chance superposition on the X-ray position. The a 
posteriori probability of this coincidence was cal- 
culated to be ~30% (based upon a 3" radius 90% 
confidence Einstein error circle). 

As part of a continuing program to probe the 
optical/UV counterparts to the luminous globu- 
lar cluster X-ray sources, we have reexamined the 
optical position in the light of new Chandra X- 
ray imaging data, and also obtained time-resolved 
HST/STIS optical spectra of the candidate coun- 
terpart to the burster X1746-370. 

2. DATA ANALYSIS AND RESULTS 
2.1. Chandra X-ray Position 

Chandra observed the field of X1746-370 with 

a short exposure of 3.2ks on 2000 May 5; we ob- 
tained the dataset from the archive. The ~30' x 30' 
field-of-view high resolution camera imager (HRC- 
I; Murray et al. 1997) was approximately centered 
on the cluster center. Data reduction was under- 
taken with routines in CIAO^. The observation was 
not affected by any periods of background fiaring, 
hence applying the standard good time intervals 

^Available at http;//asc. harvard.edu 
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and correcting for dead time, wc obtained 2903s 
of data on source. Owing to this shallow expo- 
sure, only three sources were detected at greater 
than 3(7 above background in the entire field, us- 
ing a Mexican-Hat wavelet source detection rou- 
tine (wavdetect). The brightest source (50000 
counts) was clearly the main target, whilst the 
other two were very much fainter (only 27 and 
11 net counts) and were located within 3" and 
20" of chip edges respectively, where the PSF is 
most extended and distorted. As a consequence, 
neither faint source could be reliably used to im- 
prove/confirm the nominal astrometry from the 
satellite's aspect solution, although probable opti- 
cal identifications (within ~1") to the bright (non- 
cluster) stars HD 161892 and GSCii S2220122126 
were made. For the reprocessed version of the 
dataset we used, there are no known aspect off- 
sets, hence we take the ccntroidcd position of 
X1746-370 given by wavdetect as the best X- 
ray position possible: a = 17''50™12?73±0?05, 
5 = — 37°03'06'.'8±0'.'6, where we quote our uncer- 
tainties as those of the Chandra aspect (Aldcroft 
et al. 2000). 

Wc have estimated our sensitivity limits for 
other lower luminosity cluster sources. Within the 
cluster half-mass radius of 38" (Harris 1996), but 
beyond the ^20" wings of X1746-370, a source 
with Lx ~ 10'^^ ergs^^would be just detectable, 
yielding >5 counts'^, comparable to the brightest 
quiescent LMXBs detected by Chandra in other 
clusters (see §1). Closer to the bright source the 
limits will naturally be higher due to the increase 
in effective background; this includes the entirety 
of the cluster core region. 

2.2. Refined HST Optical Position 

Although a precise optical position for the pro- 
posed counterpart to X1746-370 (star Ul) was 
published by Deutsch et al. (1998), that astrome- 
try was based on the HST Guide Star Catalog ref- 
erence frame (Lasker et al. 1990). This predates 
the improvement made possible by Hipparcos, and 
the construction of the International Celestial Ref- 
erence System (ICRS; H0g et al. 2000), the rcfcr- 

^ We assume a distance of 11.2 kpc and ATjj = 2.1xl0^^cm~^ 
(as determined from BeppoSAX spectra, Parmax et al. 
1999), and typical blackbody or thermal bremsstraiilung 
spectral models representing quiescent LMXB or cata- 
clysmic variable star emission. 



encc frame for the Chandra astrometry. Although 
the difference in a given position is typically ^ 1" , 
this is still significant when compared to the sub- 
arcsccond precision of the Chandra X-ray position. 
Hence we have redetermined the astrometry using 
the very accurate USNO-A2.0 star catalog (Monet 
et al. 1998), which makes it possible to tic an ar- 
bitrary field rather easily to the ICRS with sub- 
arcsecond precision. Once again we make use of 
the ground-based CCD image of the NGC 6441 
field kindly provided by G. Jacoby (e.g. see Ja- 
coby et al. 1997) . We select 46 bright, isolated 
stars in common between the USNO-A2.0 catalog 
(epoch 1982.0 in this field) and the ground-based 
CCD image (epoch 1995.6) and fit an astrometric 
solution to the image using IDL procedures writ- 
ten by E.W.D. and from the Astronomy User's 
Library (Landsman 1993). The residuals of the 
fit ((7=0'.' 56) imply an approximate uncertainty 
{a/y/n - 3) in the alignment to the USNO-A2.0 
frame of 0'.'09 before considering proper motion 
effects. Deutsch (1999) derived empirical imcer- 
tainties in the transfer of the USNO A-2 frame (via 
reference star matching) to ground-based images 
comprising the Second Digitized Sky Survey (DSS- 
II), which effectively includes the scatter induced 
by the differing epochs of the data and random 
proper motion effects. Hence, we conservatively 
adopt the la radial uncertainty of 0'.'35 given in 
the last row of Table 1 in Deutsch (1999). 

The next step is to transfer this solution 
to the U-band (F336W filter) images from the 
HST/WFPC2 (epoch 1994.7). These data are de- 
scribed in detail in Deutsch et al. (1998). The 
PC was centered approximately on the core, and 
hence the corresponding portion of the ground- 
based image is very severely crowded. It is there- 
fore clearly advantageous to use the data from the 
WF chips outside of the cluster core and relate 
this to PC chip on which the candidate counter- 
part star Ul is located. This is made possible by 
the STSDAS routine metric, which uses the well- 
calibrated relative positions of each of the chips 
and the geometrical distortions across the field of 
each to accurately give an RA and Dec for any star 
on any chip, based upon the astrometric solution 
present in the WF2 image header. We there- 
fore identified 52 well-isolated stars which appear 
on both the ground-based image and the 3 WF 
chips. Given the ICRS astrometry solution now 
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written into the header of the ground-based im- 
age, we are able to obtain ICRS positions for each 
of these stars and cross-compare to the nominal 
result given by metric for the WFPC2 astromct- 
ric solution. We find that shifts of 1'.'6 and O'.'OOl 
in RA and Dec are adequate to bring the WFPC2 
data onto the ICRS, with no correction to the 
orientation necessary. Again from the standard 
deviation of the residuals we estimate a radial 
uncertainty of 0'.'02 for this step, negligible com- 
pared to the other uncertainties. With the correc- 
tions in place, we measure a new position for star 
Ul in the ICRS of: a = 17''50'"12?728±0?029, 
5 = -37°03'06'.'53±0'.'35. This position differs 
appreciably (> 1") from that quoted by Deutsch 
et al. (1998), with essentially all the difference due 
to the different astrometric frames, and should su- 
percede those earlier data. 

2.3. Optical Spectroscopy 

On 1999 June 28, we obtained HST optical 
spectroscopy of star Ul. We used STIS (Woodgate 
et al. 1998) with a 0'.'2x52" aperture which was 
well-centered on the target, following a blind off- 
set from a nearby bright star. The G430L grating 
was employed, together with the STIS-CCD detec- 
tor, yielding a useful wavelength range of~3000A- 
5700A and spectral resolution « SA. Five HST 
orbits of data were taken, with five separate expo- 
sures during each, except for the first HST orbit 
for which there are only four exposures. 

The standard pipeline reductions produce a va- 
riety of data products, where all the data taken in 
a given HST orbit have been combined together. 
This combining facilitates the removal of the nu- 
merous cosmic ray events, via one-sided sigma 
clipping (together with a noise model for the de- 
tector). Unfortunately, the CCD also suffers from 
a number of defects, which appear as anomalously 
high or low spikes in the extracted spectra. These 
we identified and later removed by hand from the 
spectra. Taking the combined cosmic-ray rejected 
spectrum images from each HST ovhit, we applied 
the STSDAS xld aperture extraction and calibra- 
tion routine to produce 1-D flux and wavelength 
calibrated spectra for the target and a number of 
other reasonably isolated stars of similar bright- 
ness located on the long slit. For each star, we 
kept the width of the extraction regions the same 
for all images, but allowed the position to track 
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Fig. 1. HST /STIS optical spectra (not dcred- 
dened). Upper: from all five HST orbits of data 
- the candidate counterpart to the X-ray burster, 
star Ul (solid), with a (typically red) comparison 
star over-plotted (dashed). With approximately 
constant flux density, Ul is relatively extremely 
blue. Lower: spectra of Ul from two distinct or- 
bits when it was at its brightest and faintest ( to aid 
clarity the data have been boxcar smoothed over 
3 pixels), demonstrating the change in flux levels, 
visible even in these dispersed spectra. We also 
note the lack of appreciable changes in spectral 
slope, implying that we have obtained reasonably 
contamination free spectra. Note: the gaps in the 
spectra are regions where the data were corrupted 
by CCD defects. 

any change in the trace position. We were careful 
in setting both the source and neighboring back- 
ground regions to exclude as much as possible con- 
tributions from other cluster stars, but inevitably 
in this dense core region, many very faint stars 
will be unavoidable. Examination of these spec- 
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tra shows that: (i) Ul possesses a featureless con- 
tinuum with no strong emission hnes or notable 
absorption features (apart from Can H+K from 
imperfect background subtraction) at the modest 
signal-to-noise of these data; (ii) in comparison 
to other stars the (reddened) continuum of Ul 
has almost constant flux density, making it the 
brightest star in the blue (< 4000 A) but one of 
the faintest in the red; (iii) relative to the com- 
parison stars the flux level of the continuum of Ul 
appears to vary significantly. In figure 1 we show 
the median stacked spectra (total exposure time of 
12680s) of star Ul and a representative compari- 
son star, and two individual spectra showing the 
variability of Ul between its highest and lowest 
flux levels. As this cluster is substantially red- 
dened {Eb-v ~ 0.4), the observed flat spectrum 
implies intrinsic colors that are very blue. 

In order to compare the flux levels to the pre- 
vious WFPC2 photometry, we have estimated the 
corresponding STMAG magnitudes from our spec- 
tra: for F336W ([/) they range from 18.72±0.G2 to 
18.94 ± 0.05 and for F439W (B) from 18.76 ± 0.02 
to 18.91 ± 0.05, with uncertainties based on sys- 
tematic variations (see below). Comparing to Ta- 
ble 1 of Deutsch et al. (1998), the F336W mag- 
nitudes are consistent with the 1994 levels given 
the source's variability, but at F439W the source 
is -^0.3 mag brighter than in 1994, i.e. it ap- 
pears considerably redder. It is possible that there 
is spectral contamination present from a redder 
star (e.g. the wings of a neighbor, centered just 
off the slit, 0'.'35 away), but the lack of signifi- 
cant color-terms in the variability suggests not (see 
fig.l lower panel). Moreover, we note that the 
1995 WFPC2 data found the source with similar 
F439W brightness levels to our 1999 spectropho- 
tometry, and since no F336W data were taken in 
1995, the source could have been similarly redder 
then. 

In order to further quantify the variability ex- 
hibited by Ul, we have calculated the total net 
counts in the continua of Ul and several compar- 
ison stars for each HST orbit. We repeated the 
aperture extractions from the combined cosmic- 
ray rejected images, but this time without apply- 
ing any further calibrations, yielding count rate 
versus pixel spectra. Next we fitted the contin- 
uum shape of the source and background spec- 
tra with spline functions, and replaced points de- 



viant by more than 4(t, by the appropriate contin- 
uum fit. Lastly, we summed the counts in these 
cleaned 1-D spectra over all pixels. At this stage 
we also estimated the random errors in the net 
counts, based upon Poisson counting statistics to- 
gether with CCD read-noise contributions. These 
indicated that a high nominal accuracy can be ob- 
tained, and therefore that variations seen in the 
comparison stars (which should of course exhibit 
constant fluxes) must arise from systematic effects. 

To examine these systematic eff'ects further we 
considered higher time-resolution data. Instead of 
combining all the exposures from a given HST or- 
bit to remove cosmic rays, we used all available 
adjacent pairs (i.e. calculating a running average) 
giving 3 or 4 cleaned spectra per orbit. Once again 
systematic trends were apparent. At least two sta- 
tistically independent data points are available per 
orbit, so we use the scatter of the most divergent 
pair to conservatively estimate the uncertainty in 
the count rates we obtained from the fully com- 
bined data. The final lightcurves are plotted for 
Ul and the comparison stars in figure 2. System- 
atic long-term trends in count rate are immedi- 
ately apparent, but these are reasonably well fit 
(i.e. xt < 1) by a linear function in every case 
apart from Ul, where xt = 7. In the lower panel, 
we plot the data for Ul once more, and this time 
fit sinusoidal models to the data, one with a con- 
stant offset and another allowing for a linear trend, 
akin to that of the comparison stars. In each case, 
the period is constrained to that of the X-ray dips 
(5.73 hr), and we see that the observed optical 
modulation is indeed consistent with this periodic 
model. The fits are clearly much better than for 
the simple linear trend, with xt — 1-8 and 0.9 
for the sinusoid, without and with a linear term 
respectively. Applying an F-test to the two mod- 
els indicates that the addition of the linear term 
to the sinusoid is only marginally preferred at the 
66% confidence level. We find peak-peak ampli- 
tudes of 16% and 20%, in reasonable agreement 
with the short term ~30% changes seen previously 
in this star in the F336W photometry (Deutsch 
et al. 1998). While we certainly cannot claim to 
have detected the X-ray period in the optical data, 
we feel our data do provide evidence for this eflfect. 
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Table 1: X-ray and optical positions for X1746-370 (in the ICRS). 



Band 


Telescope 


a (J2000) 


Error 


S (J2000) 


Error 
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(la) 


X-ray 


Chandra/miC-I 


17 50 12.73 


0?05 


-37 03 06.8 


0'.'6 


Optical 


HST/WFPG2 


17 50 12.728 


0?029 


-37 03 06.53 


0'.'35 
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Fig. 2. — Lightcurvcs for the candidate coun- 
terpart star Ul and three other stars of similar 
brightness; datapoints are derived from the or- 
bitally stacked spectrum images and (systematic) 
errors estimated from higher time resolution anal- 
ysis. Upper: a linear trend has been fitted to each 
curve, but is notably a very poor fit to the data 
for star Ul. Lower: considering star Ul only, two 
models have been fitted: (i) a constant -|- sinusoid 
(dashed) with period constrained to 5.73 hr (X- 
ray dip period), (ii) the same, but with a linear 
term added (dotted). 

3. DISCUSSION 

The new Chandra X-ray and optical positions 
for X1746-370 and its proposed counterpart Ul 




Fig. 3. — Comparison of the various positions for 
X1746-370 and its optical counterpart. An F336W 

PC image from HST is shown, with the old Ein- 
stein error circle (dashed) and new Chandra circle 
(solid) overlaid (both 90% confidence). The two 
blue stars identified in the larger error circle by 
Deutsch et al. (1998) as well as their proposed UV- 
bright counterpart, Ul, are indicated. We note 
the excellent positional agreement between the re- 
duced X-ray error region and the position of Ul. 

are summarized in table 1, and illustrated in fig- 
ure 3. The positional agreement is excellent and 
well- within the (small) estimated radial uncertain- 
ties, providing strong support for the optical iden- 
tification. Following the a posteriori probability 
estimate of (Deutsch et al. 1998), but with a new 
90% confidence radius of I'.'l for the Chandra er- 
ror circle, the area enclosed is 7 times smaller 
than for Einstein, and hence only a ~4% prob- 
ability remains that we have chance alignment of 
an unrelated UV-excess cluster star with the X- 
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ray position. Moreover, our new HST/STIS spec- 
tra confirm that star Ul is unusually blue. But 
more conclusively, the temporal coverage afforded 
by our timc-rcsolvcd spectra arc sufficient to con- 
firm the optical variability of the source, which 
appears periodic and can be well fitted with a si- 
nusoid constrained to the 5.73 hr X-ray period. 

The firm identification of the optical counter- 
part does have further implications for our view of 
the nature of the source, including the issue of the 
unusual energy independent X-ray dips. First, the 
most likely origin for the sinusoidal optical mod- 
ulation on the X-ray period is the varying contri- 
bution from the bright X-ray heated face of the 
donor. Hence, its possible detection would not 
only confirm the correct identification, but also 
support the standard picture that the recurring 
dips are related to periodic obscuration on the or- 
bital period. Second, as previously noted by Par- 
mar et al. (1999), the faintness of the optical star 
implies that Lx /Lgpt ^ 1000, typical of LMXBs in 
which we directly observe the central source, and 
therefore consistent with the detection of bursts 
from this source. In particular, this suggests that 
it is unlikely that the dips are primarily due to 
obscuration of an extended accretion disc corona 
(ADC); in the classical ADC sources where only 
scattered X-rays are visible, Lx/Lopt ~ 20. Our 
optical spectral data also argue against such an 
ADC interpretation. Considering all LMXBs with 
comparable periods (and hence disc sizes), most 
show characteristic emission lines in this spectral 
region at H/?, Hen A4686 and the Bowen Ciii/Niii 
blend A4640; moreover, all three known classical 
ADC sources exhibit strong emission at one or 
more of these lines (see e.g. van Paradijs 1995). 
Hence, the lack of strong line emission in XI 746- 
370 is unlike any of these ADC systems, and is 
even unusual compared to the field LMXBs in 
general. However, wc note that our HST spectra 
of the globular cluster LMXBs in NGC 6624, and 
6712 (Deutsch 1998) are similar to X1746-370 in 
NGC 6441 insofar as they are also very blue/UV, 
but largely featureless at modest signal-to-noise 
and resolution. 

Within the context of the diverse nature of the 
cluster LMXB population, X1746-370 might at 
first glance be considered a rather average system. 
In terms of period it lies midway between the three 
with ultra-short periods (P< 1 hr) in NGC 6624, 



1851 and 6712 (see e.g. Homer ct al. 2001a, and 
references therein), and the two long period sys- 
tems in Terzan 6 and M15 (AC211) (P=12.4 hr 
and 17.1 hr respectively; In't Zand et al. 2000; Ilo- 
vaisky et al. 1993). However, it is the longest pe- 
riod confirmed burster, since the Terzan 6 LMXB 
has not been seen to burst, and we now know 
that AC211 (M15-X1) is almost certainly an ADC 
source and M15-X2, a source without a known or- 
bital period, is the probable burster there (White 
& Angelini 2001). The recent results on NGC 6652 
(Heinke et al. 2001) indicate that its burster is 
perhaps the most similar, with a longest possible 
period of 4.4 hr, though 0.92 hr also fits the avail- 
able data. If this shorter period in NGC 6652 does 
turn out to be correct, of the five bursters with de- 
termined or well-constrained periods, four would 
be double-degenerate ultra-compact systems and 
only X1746-370 would be similar to a typical galac- 
tic burster. Deutsch et al. (2000) have already 
commented on this prevalence of very exotic sys- 
tems in globular clusters. It would appear that the 
unique formation/evolution processes at work in 
globular cluster cores (see e.g. Hut et al. 1992, for 
a review) e.g. the tidal capture and exchange en- 
counter mechanisms and subsequent stellar inter- 
actions, leading to the hardening of already hard 
binaries may have led to an enhancement of the 
ultra-compact LMXBs at the expense of the wider 
systems like X1746-370. In any case, the (grow- 
ing) population of LMXBs may serve as impor- 
tant tracers of the stellar dynamics and evolution 
within globular clusters. 

We once again thank George Jacoby and col- 
laborators for providing ground-based CCD im- 
ages of NGC 6441. We are grateful to the anony- 
mous referee for useful comments. Support for 
this work was provided by NASA through grant 
NAG5-7932, as well as STScI grant GO07363.01- 
96A, and SAO grant GOO-lOllX. 
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